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We consider the general problem of estimating the inflight LISA noise power spectra and cross- 
spectra, which are needed for detecting and estimating the gravitational wave signals present in the 
LISA data. For the LISA baseline design and in the long wavelength limit, we bound the error on 
all spectrum estimators that rely on the use of the fully symmetric Sagnac combination (£). This 
procedure avoids biases in the estimation that would otherwise be introduced by the presence of 
a strong galactic background in the LISA data. We specialize our discussion to the detection and 
study of the galactic white dwarf-white dwarf binary stochastic signal. 
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I. INTRODUCTION 



The LISA mission is a collaboration between ESA and NASA to build and operate a space-based laser interferometer 
aimed at detecting and studying gravitational waves (GW) of astrophysical origin. The current schedule places the 
launch around theyear 2011, and the 3-10 years observation period is expected to produce a number of observations of 
galactic binaries [l|, of massive black hole mergers and captures, and possibly of a cosmological stochastic background 
I | • Early studies of the data analysis techniques that will be required to observe these signals are important since 
^ ■ they couple our present knowledge about the sources LISA will be able to observe to the science requirements guiding 
the design of the LISA mission itself. It is in this context that we present here a study of the accuracy that can be 
expected for the characterization of the noise in the LISA detector once it has been placed in its heliocentric orbit. 
This is needed not only for being able to assess the LISA sensitivity in a region of the frequency band where the 
stochastic signal from coalescing galactic binaries would prevent it, but most importantly for being able to detect 
weak gravitational wave signals with confidence. Theoretical modeling and pre-launch performance measurements 
of flight-unit instruments will certainly be performed. However these might not be sufficient to achieve the level of 
• accuracy required in the analysis of the data collected during the inflight operation of LISA. 
q-i Specifically, we calculate in the low frequency regime the optimal approximations for the spectra of the noise affecting 
I . the interferometric combinations sensitive to gravitational waves, by relying only on all the possible cross-spectra 
between the interferometric combinations that are sensitive to gravitational radiation and the symmetrized-Sagnac 
combination, £. The rational behind this is that £ couples very weakly to gravitational radiation while it is affected 
by the same instrumental noises as the other combinations [Tj. These approximations give a lower bound on the 
error that can be achieved by estimators which make no prior assumptions about the characteristics of the various 
noise sources. Such estimator could be used for carrying out the full analysis, or for validating the noise models one 
could adopt for implementing more sophisticated and efficient data analysis techniques. An outline of the paper is 
given here. 

Section [n] presents a review of the LISA baseline design, which is then followed by a discussion of the noise 
characterization problem in section llTll In scction lTvl we show that closure for the solution to the noise characterization 
problem does not exist, and that only approximate solutions can be implemented. In section fVl we describe the noise 
model used in our calculations, which accounts for possible cross-correlations between pairs of various noise sources. 
Although these terms have been neglected in the literature, they could contribute to the overall noise budget of the 
mission, and for this reason they are included in our analysis. We then present in section IVll our calculation of the 
lower bound on the estimated noise spectra in the LISA interferometric combinations that are sensitive to gravitational 
radiation. In section IVlII we finally provide our comments on the implications of these bounds for the astrophysical 
reach of the LISA mission, and on possible implementations of noise estimation algorithms that would approach the 
error bounds we estimate. 
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FIG. 1: Schematic description of the optical benches on spacecraft 1 in the LISA baseline design. Variables on the right hand 
bench are distinguished from corresponding variable from the left hand bench by asterisks. The measured Doppler signals are 
2/3i and 1/21 • nt denote normal vectors along the triangle arms, and V\ denotes the random velocity vector of proof mass 1. 

II. LISA BASELINE DESIGN AND NOTATION 

The LISA mission will consist of three spacecraft flying in an approximately equilateral triangle formation. The 
length of the arm opposite to spacecraft i (i = 1,2,3) is labeled Li. The three nominal arm lengths will be equal 
to 16.7 s (c = 1), and they will differ from each other by at most a few percents during the entire duration of the 
mission. In contexts where the differences between the length of the arms can be disregarded, we will simply refer to 
the nominal arm length L, without indices. As it is illustrated in figure ^ each spacecraft will contain two optical 
benches. Each bench will be designed to bounce a laser beam between each of its proof mass and another proof 
mass on a distant spacecraft. The proof masses are mechanically isolated from their optical bench, and their relative 
position is measured with electrostatic sensors. These generate inputs for the drag-free control systems that allow 
the spacecraft to isolate the proof masses from external disturbances in such a way to maintain them (to a very 
good approximation) in a free-falling configuration in the directions to the other two spacecraft. In addition, laser 
beams will be bounced off the back of the proof masses and will be exchanged between the optical benches on a given 
spacecraft in order to sense the motion of the optical benches relative to the proof masses. 

With this baseline design of the LISA spacecraft, the frequency noise of the lasers will dominate the secondary 
noise sources (such as acceleration noise of the proof masses or optical path length distortions) by approximately 
seven orders of magnitude in the frequency band of interest (10~ 4 — 1 Hz). Time delay interferometry (TDI) has 
been proposed 0, HQ as a robust way for suppressing the frequency fluctuations from the lasers and the mechanical 
vibrations of the optical benches to levels smaller than that identified by the secondary noise sources. This is done by 
properly time-shifting and linearly combining the measurements of the phase differences between the coherent laser 
beams exchanged between the three pairs of spacecraft and the three pairs of on-board optical benches p| . In reference 
it has been shown that the entire space of the TDI combinations can be generated by four combinations, called 
a, j3, 7 and £, which are the Sagnac Interferometric combinations LISA will be able to synthesize. Since (a, (3, 7, £) can 
be regarded as a basis for the entire space of the LISA TDI combinations, in what follows we will limit our analysis 
to these four interferometric data. 

It should be noted that the modifications of the Sagnac observable (a,/3, 7, C) required in order to account for the 
rotation of the LISA array and variation in time of its arm lengths Q can be disregarded in the analysis presented 
here. Although it will be critical to account for these effects in order to suppress the laser phase fluctuations below 
the secondary noise sources, the modifications introduced to the gravitational wave and secondary noise responses 
entering into the TDI combinations a, (3, 7 and £ (derived under the assumption of a stationary LISA configuration) 
will be negligibly small. 

It is anticipated that the dominant noises affecting all the TDI combinations will be introduced by the proof 



3 



masses and optical path noises (shot noise and beam pointing fluctuations). In the notation of Fig. ^ the frequency 
fluctuations generated by the mechanical vibrations of the proof masses are given by || 

unroof mass ~ / „. , „.* „.* \ 

C = "1 • (^2,2 - V 2 ,1S + V 3 3 - W 3)2 l) 

+fl 2 ■ (v 3 ,3 ~ V 3 .21 + Vl ± - V* t23 ) 

+n 3 ■ (vi,i - vi,23 + v* 2 2 - vl 13 ), (1) 

where bold characters are either vectors or matrices, Vij = Vi(t — Lj), and Vijk — v i(t — Lj — L^)- We shall use 
below the definitions v± — n 3 ■ v 3 , etc. Denoting by rnj the optical path noise affecting the Doppler measurement yij, 
the optical path noise contribution to ( is given by the following expression p| 

^ptical path = ns2 2 _ ^ + ^ _ nsi i + n2i i _ ^ (2) 

If we denote by £ GW the contribution to the ( observable from a gravitational wave yj Eq. (60)] , it follows that ( 
can be written in the following form 

C = C GW + C ins , (3) 

where the instrumental noise part is given instead by 

^ins ^proof mass _j_ ^-optical path /^\ 

In what follows we will be interested in calculating the spectrum of the £ combination 

S (( = E[CC], (5) 

where £ is the Fourier transform of £, and E denotes the expectation value over many realizations of the instrumental 
noise. In the particular case of a stochastic gravitational wave signal, the expectation value will be performed over 
the realizations of the signal too, although it always will be assumed to be independent of (and therefore uncorrelated 
with) the noise. Under these assumptions we can write 

SqQ = S^GW^GW + S^ins^ins . (6) 

In what follows, for any pair of interferometric observables (x,y), S xy will refer to the following function of the Fourier 
frequency / 

S xy = S* yx = E[xy*}- (7) 

For simplicity of notation we will write S®y for S x Gw y Gw and S x n y for S^ms,,™. As implied by Eas. HH2|) . the spectrum 
can be expanded in terms of proof mass and optical path noises as follows 

T x 3 2 3 



S™ = 4 sin 2 ( - Y^[Sv iVi + S vjvj + 2Re S ViVJ ] +Y, + Sn ^J 

' i—l i—1 j=i-\-l 

— 2Re [S ni2 , ll3 + S n2in23 + S n31 n 32 ], (8) 



where Re x denotes the real part of x (Im x is the imaginary part), and where for the sake of generality we have 
included terms accounting for correlations between proof masses on-board the same spacecraft and between optical 
path noises along the same arm (see section [V] for the physical motivations behind the inclusion of these correlations 
in our analysis). 

III. IN-FLIGHT NOISE CHARACTERIZATION 

As it is well known, the optimal statistics to detect a signal s in noisy data is the matched or Wiener filter, at least 
when the noise distribution can be well approximated as Gaussian. For the combination a, for instance, this statistics 
is 

df^, (9) 
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where s a is the signal s filtered by the transfer function of gravitational wave to the combination a. As shown in 9], 
the three TDI Sagnac observables a, (3, and 7 can be combined to form the new data set A, E, and T. The advantage 
of these three combinations is in their noises being statistically uncorrelated for every frequency /, 



E[AE*} = 0, 

and similarly for all possible pairs taken among A, E, and T. 
diagonalizing the correlation matrix 



(10) 

The uncorrelated combinations are obtained by 



Solo. 
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(11) 



leading to an optimal detection statistics for the signal s that is given by 

J dfxC- 1 ^, (12) 

where x — [a, (3, 7] and x — [sa, S/3, s 7 ], for s a , sp, and s 7 the signal s filtered by the transfer functions of gravitational 
waves to the combinations a, (3, and 7 respectively. Eq. is only a special case of Eq. I|12f> . with C, x and % 
restricted to the one-dimensional case. 

The optimal signal-to-noise ratio po that can be achieved by Eq. Ijl2|l is 9] 



Pl = J dfxC- 1 ^ 



(13) 



Any variation of the detection statistics (up to a scaling factor) will give a smaller signal-to-noise ratio, so Eq. i|12|) 
is the optimal detection statistics. It is also true that the detection statistics of Eq. 1)12(1 maximizes the probability 
of detecting the signal for a fixed value of the false alarm probability. 

Since the correlation matrix can only be known with some finite accuracy, we can write it in the following way 



C = Co + Ci, 



(14) 



where Co is the true correlation matrix, and C\ is the error in its determination. By definition, the matrices C, 
Co and Cl must all be hermitian and positive definite. When the correlation matrix is measured with some error, 
the detection statistics of Eq. I|12l) constructed with the measured matrix is no longer optimal with respect to the 
true correlation matrix, and the SNR can be expected to be reduced. We present in Appendix lAl a derivation of the 
SNR loss due to the presence of the error C\. To give an example of the magnitude of this effect, we can assume 
the estimate of the correlation matrix to have a given relative error, (COy = e(Co)jj. Considering a fairly general 
signal [cf. Ea. (|A6J) ]. the relative SNR loss is then ~ e/2. These considerations indicate that in order to limit the SNR 
reduction to a negligible level (say ~ 0.1%), the relative error on the measured spectra and cross-spectra for the three 
TDI combinations a, (3, and 7 should be of order 10% at most. 

Errors in the determination of the spectra of the TDI combinations will also have an effect on the measurement 
of the parameters of the observed gravitational wave signals. Let J(A) be the Fisher information matrix for the 
estimation of the parameters A using data x with distribution function p(x; A): 



/(A) = -E 



d 2 \ogp(x- A) 
d\ 2 



(15) 



The Cramer- Rao bound gives a limit on the error cov A of an unbiased estimator \{x): 

cov A = E[(X(x) - A) t (A(a;) - A)] > / _1 (A). (16) 

As usual, the > sign in Ea. HRj|) means that the left hand side of the inequality minus its right hand side is positive 
semi-definite. An unbiased estimator can achieve the Cramer-Rao bound if and only if it has the form 



\(x) = A 



<91ogp(a;; A) 



1 T 



I" l (A) 



(17) 



Specializing to the case relevant here, where x = [a, (3, 7] and where x(A) = [s a (A), s / g(A), s 7 (A)], assuming Gaussian 
noise, and going to the Fourier domain, one can write Eq. l|17|l into the following form after some simple algebra 



\(x) = A + Re 



(18) 
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where 

This estimator achieves the Cramer- Rao bound only if the right-hand side of Ea. i|18|) is independent of A, as the 
notation suggests; if this is not the case, no unbiased estimator can achieve the Cramer-Rao bound. For the analysis 
of data from LISA, one will use the estimator in Ea. (|18fl . with some estimated value of C that contains an error, 
as in Ea. (|14f> . This estimator will necessarily have larger errors than the Cramer-Rao bound. Because the second 
term of the right-hand side of Ea. ((T%|l is linear in x — x, X(x) will be unbiased, independently of the size of the 
error C\ . It is shown in Appendix ^ that the error term in the covariance of the estimator constructed with the 
perturbed correlation matrix will vanish to first order in the perturbation C\, so that, like the SNR, the variance of 
the estimated quantities will be at most second order in the measurement errors of the correlation matrix. 

At high frequencies (/ > 1 mHz, |IcjL the gravitational wave sources are expected to be discrete in time or in 
frequency, so that it should be possible to identify portions of the time series that are free of signal, and can therefore be 
used to estimate the spectra with an excellent accuracy (which is limited only by the integration time, and indirectly 
by the stationarity of the noise). At lower frequencies, on the other hand, gravitational waves from the galactic 
binaries will not be distinguishable, and will form a stochastic signal above the instrumental noise. Strong signals 
that are above this confusion limit might be detected and analyzed with respect to the noise from the binaries, which 
can be estimated with high accuracy by integrating the data for a long enough time. The most direct application of 
the characterization of the noise at low frequencies is thus the measurement of the stochastic background itself. In 
other words, it is highly desirable to have techniques that can effectively "turn off" the gravitational signal, in order 
to properly characterize the instrumental noise. At higher frequencies, these techniques should not be crucial for the 
detection of weak signals, but should be very useful in understanding and validating the detector performance. To 
give a specific example, it might be useful at some point to verify that a given spectral line is of astrophysical origin 
and is not instrumental, without having to rely on astrophysical assumptions (such as the Doppler modulation of the 
signal) . 

In the simplest search, a stochastic background will be detected by comparing the measured spectra S aa , Spp, 
and S 77 to the estimate of the instrumental noise. Differences between the measured and the estimated spectra will 
be interpreted as a stochastic background. It therefore goes without saying that the low frequency portion of the 
instrumental noise spectra in those TDI combinations must be estimated as well as possible. 

As it has been recognized in |ll|. the TDI combination ( (i.e., the fully symmetric Sagnac combination) is much 
less sensitive to gravitational waves than the other Sagnac combinations (a, (3, 7) in the low part of the frequency 
band accessible to LISA (/ < 30 mHz), although it is affected by the same instrumental noise sources. Consequently, 
it can be used as a "gravitational waves shield" to estimate the instrumental noise. This approach was introduced 
by ^3 to bound the power in the stochastic background. Since the GW signal is additive in the noise of the TDI 
combinations, 

Sxx = Sx s x + S xx , (20) 

where Sxx is the measured spectrum of the Unequal- Arm Michelson combination A, S^x ^ s the spectrum of the 
instrumental noise in X, and S'xx i s the spectrum of a stochastic background. It is shown in [TIJ that the following 
inequality holds 

S X x > s xx - Sxx, ( 21 ) 

where the upper bound on the instrumental noise spectrum is, in the low-frequency region of the LISA band, 

S^ s x = 16% - l28ir 2 f 2 L 2 V - 32(3 - 4ir 2 f 2 L 2 )N, (22) 

with V being a lower bound on the spectrum of the acceleration noise of all proof masses, and N a lower bound on 
the spectrum of the optical noise of all optical paths. 

A somewhat similar approach is used in |l2j | to estimate the average of three TDI combinations, S = (Sxx + Syy + 
Szz)/3. They estimate S lns by multiplying by a judiciously chosen transfer function, and subtract the result from 
the measured S to get S GW . This "averaged" estimator is quite different from what we will consider below, and will 
not be discussed further. 

IV. MODEL INDEPENDENT APPROACH 

At low frequencies, the real function and the three complex functions S a £, Sp^, and 5* 7 ^, give seven measurements 
that are insensitive to gravitational waves, in the sense that SgQ ~ <Sg" s , 9 = 01,(3, 7, C- Since the four TDI generators 
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fulfill the following identity 



j(Li+L 2 +L 3 )j£ _ 



a iuj(L2+L 3 ) 



e iuL 2 _ e iu(Li+L 3 ) 



(3- 



JuL 3 _ iu)(L 1 +L 2 ) 



(23) 



the seven functions introduced above are not all independent. In particular, we could try to estimate by a linear 
combination of the cross-spectra S a (, Sp^, iS 7 f, and S 1 ^ as 



Using the identity Eq. Ij23|l to expand the right-hand side of Eq. (|24J) and collecting terms in gives 
[coswii — cosuj(L 2 + L 3 )]a + [1 — coso;(i 1 — L 2 — L 3 )]z = 1 — cosa>(Lx + L 2 + £3). 



(24) 



(25) 



Further requiring all the other terms in the right-hand side of Ea. (|24() to be zero gives eight other equations. However, 
the only solution satisfying those eight equations is 



(a, 6, c, z) = (cos lo(Li — L 2 ) — cosojL 3 , i — coslo(Li — L 2 + £3), 
[coslu{L 1 + L 2 + £3) - 1] sin|(Li - L 2 + L 3 



cos u)L\ + cosuj{L 2 — L 3 ), 



sin|(Li +L 2 + L 3 ) 



(26) 



and this solution does not satisfy Ea. (|25|l . This implies that cannot be reconstructed from measurements of 
spectra that are insensitive to gravitational waves. The argument can be extended to other spectra with f3 and 7. ft 
is therefore necessary to consider combinations of the insensitive spectra that are only approximations to the spectra 
involving a, (3, and 7. 



NOISE MODEL 



In what follows we first define a realistic noise model for LISA. The best combination of the spectra and cross- 
spectra involving £ for this model will be derived in the next section, where it will be used to set a bound on the 
accuracy that can be achieved in characterizing the LISA noise spectra 

The principal source of noise at low frequencies is the acceleration noise of the proof masses along their sensitive 
axes. A number of effects lead to acceleration noises which sum up to a noise budget (per proof mass) of 3 x 10 -15 
m s~ 2 Hz -1 / 2 at 0.1 mHz 13]. The principal sources of acceleration noise might be due to thermal distortion 
of the spacecraft, cross-talk between the other degrees of freedom and the sensitive axis, gravity gradient noise due 
to spacecraft displacements, fluctuating spacecraft magnetic fields, thermal noise, back-action from sensing, electric 
forces from fluctuating charge, residual gas, magnetic damping, temperature variations, etc. It should be noted that 
the first four noise sources might affect both proof masses on a given spacecraft, so that the acceleration noises of the 
two masses will be correlated to some extent. This might be important for the correct estimation of the noise level 
that will be achievable with TDI In this paper, we will assume that the acceleration noises, converted in terms 
of relative frequency fluctuations, are given by 



)ap 
1 



^ViVi Sqj* d * 2.5 X 10 



-48 



1 Hz 



Hz" 



(27) 
(28) 



where i, j — 1, 2, 3, 5ij is Kronecker's delta, and £ is an arbitrary complex number with norm |£| < 1. The correlated 
noises should couple to the proof masses either gravitationally or magnetically, so the phase angle of £ will only depend 
on the spatial distribution of the perturbing field with respect to the sensitive axes of the two proof masses. 

At frequencies above ~ 5 mHz, the optical path noises dominate the total sensing noise. The optical path noise 
budget [13| includes contribution from the shot noise affecting the phase measurement at the photodiodes, as well as 
phase fluctuations resulting from scattered light effects, laser beam pointing instabilities, and path length variations 
on the benches and within the telescopes. The latter two effects will couple the optical path noise along various 
paths: pointing instabilities will correlate the noise along the paths related to benches i and i*, while the path length 
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variations will mostly introduce correlations at times offset by the one-way-light-time L between the measurements 
from beams going in opposite direction along a given arm. In what follow we will make the assumption that the latter 
effect is dominating the cross-correlation of the optical path noise, and write the following expressions for the optical 
path noise spectra and cross-spectra 

5 nyno . = 1.8 x 10- 37 (jLSj Hz" 1 (29) 
S nijnik = -2 X cos(27r fL)S ni . nij (30) 

where i, j, k = 1, 2, 3, i ^ j, i =/= k, j k, and x is a positive real number. The cosine term in Ea. H3U|) comes from the 
time delay L for the propagation of phase noise along one of the LISA arms. 

In what follows we will use two noise models. In the first case we will be optimistic and assume all cross-correlation 
terms to be negligible, i.e. we will pick £ = \ = 0- I n the second, pessimistic case, significant correlations will be 
assumed to be present. For the acceleration noise, we take |£| = 0.25, corresponding to the fact that roughly 50% of 
the noise budget (in strain) corresponds to noise sources that can induce strong correlations between the two proof 
masses on-board one spacecraft. Since there is no specific reason for preferring a particular value of the phase for £ 
over any other, we will treat it as a free parameter in our analysis. We will also pick \ — 0.05, based on the fact 
that roughly 15% of the optical-path noise budget (in strain) is attributed to path length variations on the optical 
benches and the telescopes. It should obviously be noted that these noise models will probably not give a realistic 
description of the noise that will affect the LISA measurements. They are only chosen to provide some guidance in 
the noise characterization problem. 



VI. ERROR BOUND ON THE ESTIMATED SPECTRA 



In this section we calculate the smallest possible errors between the estimated and real instrumental spectra that 
can be achieved by using only the GW-free measurements SgQ, where 9 = a, /3, 7, £. These spectra are the only 
ones that are free of laser frequency noise and of GW signals, and therefore form a natural basis for building noise 
estimators that do not rely on unverifiable assumptions about the properties of the noise sources. Our error bounds 
will obviously depend on the noise model used, and so will the estimators that achieve these bounds that we will 
present. These estimators will therefore not be valid estimators for the characterization of the LISA noise; model 
independent estimators that approach our error bounds will be presented in a follow-up communication. 

Each TDI combination will show a different transfer function for the noise contributions identified in the previous 
section. In particular, every spectrum will be a linear combination of the acceleration and optical noise spectra, so 
that any instrumental noise spectrum S 1 ™ (//, v are any TDI combinations) can be written as 

= T n» ' &n, (31) 
where T^ v is a complex vector of frequency dependent coefficients, and Sn is the real vector of the 21 noise spectra 

Sn \$vivi ; S v * v * , Re S VlV * , Im S VlV * , 

^v 2 V2 5 v 2 ' ^9 5 ^^2 V 2 ) 

^3 5 ^3 1 ^^3 V 3 ) IlTL Sy 3 y * , 

*S'n2in21 1 , ^Vl23«23 1 ^"21123 ' 

Sn 31 n 3 n S n32n32 , S n3in32 ). (32) 

The Sec, spectra only have terms of the form S ViVi + Re S ViV * , S v * v * + Re S ViV * , or Im S ViV '. This is because 
( has the same transfer function to Vi and to v* , so any spectrum Sg^ will only have terms of the form S ViVi + 
S*. vf and S v * v * + S*. vf . Consequently, only three of the four quantities S ViVi , S v ? v *,He S ViV * and Im S ViV * can be 
measured simultaneously. A similar argument applies for the optical noises. In that case, a Sg^ spectrum measures the 
combinations S nijnij - S*. kn and S ntkrHk - S niknij , so one can only measure S rizj n zj - Re S nikrizj , S rHk n ik - Re S nikUi] , 
and Im S nikni simultaneously. Inequalities in the lengths of the arms may break these symmetries. 

The coefficient vector T^q is known from the definition of the TDI combinations, and it just consists of geometrical 
transfer functions. Clearly, there is not enough information to measure every one of the 21 terms in Sn- This means 
that an approximated solution has to be used. Supposing we want to measure the instrumental spectrum S 1 ™, we 
can measure S a £, S 1 ^, S 1 ^, and S^, and then approximate T^„ with the following expression 

%v = aT a( + bT fK + cT 7C + zT cc . (33) 
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The estimated spectrum S 1 ^ is then given by 

= Tfu, ■ S N (34) 
= aS a £ + bSp£ + cSy£ + zS^ . (35) 

Performing the approximation in Ea. H33(l requires the definition of a proper error metric. In this paper, we choose a 
model dependent metric in order to bound the achievable error on the estimated spectrum. For complex T^„ vectors, 
we minimize the error on Re[T Ml/ ] and ImpT^] independently. Specifically, we minimize the error ||Re[T^ w ] — RefT^]!! 2 
and HlmpZ^] — Im[T MI/ ] || 2 , where Re[T M „] and ImpT^,,] are restricted to the space spanned by the vectors Reflj^] and 
ImjT^] respectively, and where 

\\xf = (x,x) (36) 

for the scalar product 

21 

(x,y)=5^a; < (S JV )?y i . (37) 
»=i 

This error norm penalizes error on individual noise components proportionally to their contribution to the spectrum 
S 1 ^. The approximation casted in this form is a linear least-squares estimation problem, which can be solved easily 
using standard techniques. In particular, the vectors T a ^, T/jf , T 7 £, and are not linearly independent, and therefore 
the solution T M „ is not unique. We pick the solution with the smallest Euclidean norm, as in [l6L Eq. (2.6.7)]. 

In order to explore the influence of correlations on the measurement of the spectra, we consider the two choices of 
£ and x defined in section For the case of no correlations, we show in figure [21 the overlaps 

/ cms C \ 

[Saa ' Sec] = 7m¥mw- (38) 

Below ~ 2 mHz, none of the Sg^ have very good overlap with S™, but above that frequency, is very close to 
S^(. We can thus expect a good least square fit for / > 3 mHz, which is what Fig. OUshows. The relative least square 
error is 



(39) 



and the relative error e on the spectrum is 



cms cms 

Z2£L Z°Z°L (A(]\ 

Cins ' \ J 

u aa 



The spectrum error peaks at -24.5% at 0.76 mHz, although it is small in some low-frequency bands. The other 
non-trivial spectrum for £ = x = is Re S^S. All cross-spectra are real, and spectra involving /3 and 7 are identical 
t° ^aa 01 ^o/f- The over l a P s are shown in FigQ] while the relative errors are shown in Fig[SJ Again, the overlaps 
are good for / > 1 mHz, and the least square and spectrum errors are small. The spectrum S 1 ^, however, has a 
zero around 0.22 mHz, which cannot be fitted by the Sg^ measurements. This leads to the large relative error on the 
spectrum below 1 mHz. 

For the pessimistic case with |£| = 0.25 and x = 0.05, the characteristics of the approximation varies with the 
arguments of £. As argued in section a zero argument for x is probably realistic. To get a specific example, we 
also fix the argument of £ to be zero, and plot the results of the least square fit for S 1 ™ in Fig. It can be seen that 
at high frequencies the error is increased significantly by the small amount of correlation in the optical path noise. At 
low frequencies, correlations do not appear to significantly increase the error of the fit. Fig. shows how the error 
varies with the phase angles of £, thus revealing the importance of the acceleration noise correlations. Depending on 
the argument of £, the estimation relative error can vary by as much as 10% at low frequencies. 

In general, some cross-spectra might be complex, but for the noise model of section [3 they will all be real, if the 
phase of £ is zero. We present in Fig. [S] the results of the least square fit for Re S^i, for the case £ = 0.25 and 
X = 0.05. As it can clearly be seen, the correlations exacerbate the problem of fitting the zero near 0.22 mHz. 
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f(Hz) 

FIG. 2: Overlaps between the instrumental noise S'^a and the measured spectra S^. By symmetry, the overlap [SJJS, Re(# 7 f )] 
is identical to the overlap [S" 1 ^, Re(S/3^)] (third panel). Overlaps with the imaginary parts of the complex cross-spectra are all 
zero. 
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FIG. 3: Results of the least square fit of S™- Top left: the weights a (continuous line), b (dashed line), c (dashed line), and z 
(dotted line). Top right: the normalized least square error. Bottom left: S™ (dotted line) and S™a (continuous line). Bottom 
right: the normalized error on the estimated spectrum. 

VII. DISCUSSION 

We have presented an analysis of the prospects for a precise characterization of the noise of the LISA detector once 
it begins to operate as a gravitational wave detector. We have derived a set of bounds on the error of the spectra 
and cross-spectra of the TDI combinations that will be necessary for the analysis of the data. The most important 
feature of these bounds is that they are derived almost without any assumptions about the form of the individual 
noise sources. The only exception is that we assume certain cross-correlations to be negligible in order to keep a 
manageable number of terms; there are, however, no fundamental limitations to including additional terms in future 
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FIG. 4: Overlaps between the instrumental noise Re S^S and the measured spectra S^. By symmetry, the overlap 
[Re S^Jf, Re Spt;] is identical to the overlap [Re i§Sf, Re S^(] (third panel). 
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FIG. 5: Results of the least square fit of Re S 1 ^. Top left: the weights a (continuous line), b (continuous line), c (dashed line), 
and z (dotted line). Top right: the normalized least square error. Bottom left: Re S 1 ^ (dotted line) and Re S a p (continuous 
line). The dips in these curves are unresolved zeros. Bottom right: the normalized error on the estimated spectrum. 

analyses. The bounds that we obtained simply reflect how well a given spectrum can be approximated by a linear 
combination of spectra that are insensitive to gravitational waves, and therefore only sample the instrumental noise. 

Our error bounds are different from those given in [13 and [13, which assumed prc-launch measurements or 
estimates of the acceleration and optical path noises to be available in constructing their estimators of the spectra. 
In contrast, we have calculated the optimal error bound that can be achieved by an estimator that does not use 
prior assumptions about the spectra of the various noises. These bounds will be most useful for the development 
of model-independent estimators, which will be the subject of a follow-up paper. It is important to emphasize that 
model-independent estimators will be critical for the unambiguous detection and study of GW sources with LISA, 
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FIG. 6: Same as FigEJfor £ = 0.25 and \ = 0-05. 
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FIG. 7: Normalized error on the estimated spectrum, for different values of the arguments of £. The continuous line is for 
arg£ = 0, and the dashed line is for arg£ = tt. These two curves bound the variation of the error at all frequencies above 0.39 
mHz. Below 0.39 mHz, a larger error (dotted line) is possible for a complex £. 



and possibly for the validation of the instrument and of other model-dependent estimators. 

To give a specific example, we plot in Fig. [§| our best-fitted spectrum and the estimated spectrum from for 
the S^x spectrum, assuming no correlations (x = £ = 0). As shown in Eq. I|22|l . the spectrum estimator from 
requires a knowledge of how small the acceleration and optical path noises can be. Without reliable prior information, 
the lower limit can only be assumed to be zero, and the spectrum estimate is S^x = 165^^, which is considerably 
worse than the least-square spectrum. One has to be able to guarantee that the noise is at least larger than 90% of 
the design value so that the spectrum estimate performs as well as the optimal least square estimator. 

If the real stochastic background is as strong as the one plotted in Fig. it will be detectable unambiguously by 
just comparing the measured spectrum Sxx to sixteen times the measured spectrum S^. However, the measurement 
of the excess power at each frequency of Sxx will carry a systematic error from the imperfect characterization of 



12 




FIG. 8: Same as Fig. [S]for £ = 0.25 and x = 0.05. 



the noise. In the model that we considered above, this error would bring down the largest frequency at which the 
background can be reliably estimated from ~ 2 mHz to ~ 1 mHz. 

To conclude, our bounds show that the spectra relevant to LISA's data analysis can only be measured to relative 
errors of the order of 10% if the estimation is restricted to using only the model- independent, signal-independent 
quantities Sq^, for 9 = a, 0, 7, £. It should be possible to unambiguously detect large stochastic backgrounds, by 
comparing Sxx to 165^, for instance. Uncertainties in the characterization of LISA's noise will principally introduce 
errors in the estimation of the power spectrum of the stochastic background. 
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APPENDIX A: SNR LOSS AND DEGRADATION IN PARAMETER ACCURACY 

We present in this appendix a derivation of the formulae describing the degradation of the SNR and of the error on 
the estimated parameters of a GW signal due to errors in the estimation of the spectra and cross-spectra describing 
the noise in the LISA detector. For matched filtering, the detection statistics is given by 

A = y dfxC- 1 ^. (Al) 

The SNR p for this statistics can easily be shown to be 



~ f 

pZ JdfxC^CoC-^' (A2) 

where the correlation matrix is assumed to contain an error as described by Eq. I|14|) . When the correlation matrix 

is perfectly known (Ci = 0), Ea. (|A2|) reduces to Ea. l|13|) . The Co matrix in the denominator of Eq. l|A2(l comes from 
the expectation of |A| 2 when only the instrumental noise is present. 
If C\ <C Co, the Neumann expansion can be used to write 

C- 1 = Co 1 - C^CiCo 1 + C^dC^dC^ 1 + 0(Cf). (A3) 
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FIG. 9: Top: a comparison of our least square bound (thin continuous line) with the bound of [TH Eq. (12)] (dash-dotted 
line for a lower bound on V and N (S° and S 1 in the notation of [TH ]') of 90% of their design values, and the thin dashed 
line for a lower bound of 50% of their design values), for the sensitivity of Sxx (dotted line) with £ = \ = 0, for 1 year of 
observation and a signal to noise ratio of 5 (the sensitivity is averaged over sky position and polarization angle, as in pj). 
The thick continuous line shows the sensitivity for 16^^, the limit of the method of [Hj when V and N are allowed to be 
vanishingly small. The thick dashed line shows the confusion noise from binaries calculated by fT^] (dashed line in their Fig. 
2). Bottom: the normalized systematic error on the sensitivity of X (continuous line, left axis), computed from the least square 
bound minus the sensitivity of Sxx > divided by the strain in the stochastic background, and the ratio of the sensitivity of Sxx 
to the strain in the stochastic background (dotted line, right axis). 



Expanding Ea. i|A2|) in the small matrix Ci, straightforward but lengthy algebra shows that the first order terms 
vanish, and that 



2 2 

P = Po 



Po J Po \J 



(A4) 



where po is the SNR achievable when C\ = 0. 

If C is known up to a scale factor for all frequencies where the signal is present, for instance, then Ea. (|A4|l shows 
as expected that the SNR is unaffected. If, however, the error on C is frequency dependent over the frequency range 
where the signal contributes significantly to the SNR, the filter used to construct A will be far from optimal, and 
significant degradation in the SNR will occur. Supposing for simplicity that the signal contributes uniformly to the 
SNR over a bandwidth B, and supposing that the error on the correlation matrix is C\ = eCo over a bandwidth b 
included in the signal frequency range, then 



pI 



i 



B 2 



(A5) 



The SNR losses will be the smallest when the error on C is uniform over most of the bandwidth of the signal, as it 
is the case for narrow-band signals, and obviously when these errors are small where the signal is large compared to 
the instrumental noise. The worst case scenario is a signal with b/B = 1/2, in which case the fractional SNR loss is 

(A6) 

As similar perturbation analysis can be performed for the parameter estimation problem, by expanding Eas. H18|) 
and (|19fl to first order in Ci, and then computing cov A as defined in Ea. (|16fl . The inverse of the perturbed Fisher 
matrix given in Ea. is given by 




lo 1 



Cn 1 C\Cr ) 



0(C?), 



(A7) 
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where Iq is the achievable covariance when the correlation matrix is known exactly (Ci = 0). Using this and the 
first order expansion of the inverse of the correlation matrix into Eci. (|18|) gives four terms that are first order in C\. 
These four terms can be simplified using the following identities 

E[(x - xY(x - x)} = Co (A9) 
E[(x - X f(x - x)\ = 0. (A10) 

The result is that the terms that are first order in C\ all cancel, as they did for the SNR. The increase in the 
covariance of the estimated quantities is therefore at most second order in the errors on the measurement of the 
correlation matrix. 



Rc 



Re 



~ dx 



C- L (x-x) 



_At 



(A8) 
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